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Modern redshift surveys such as the 2 degree field Galaxy Redshift Survey (2dFGRS) and the 
Sloan Digital Sky Survey (SDSS) reveal the fully 3 dimensional distribution of a million or so 
galaxies over a large cosmological volume. Visually galaxies appear to be distributed along sheet- 
' like and/or filamentary superclusters. The CfA Great Wall, Southern Great Wall and the recently 

discovered SDSS Great Wall are very spectacular superclusters. Clearly theoretical predictions 
for galaxy clustering must be tested against these rich datasets. This can be achieved by means 
of the Minkowski Functionals (MFs). A MF-based approach provides an excellent description of 
superclusters and voids and allows one to quantify the properties of the cosmic web. In this review 
we give a summary of the progress made in this direction. After reviewing the status of observations 
and of numerical simulations, we comment upon the nature of bias which serves as a link between 
theoretical predictions and observations. We also summarise the methods developed for efficient 
. numerical estimation of MFs for cosmological datasets and list several important results obtained 

■ using these methods. Specifically, we stress the discriminatory power of MFs and of the derived 
ly^ | morphological statistics, the Shapefinders. Shapefinders are an excellent tool with which to study 

the shapes and sizes of superclusters and voids. We also discuss some of the important effects 
of scale-dependent bias which are brought out by a MF-based study of the mock catalogues of 
galaxies. Such effects, we note, should be accounted for before comparing theoretical models with 
observations. 
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INTRODUCTION 



During the past decade, CMB experiments such as WMAP type-la supernova searches and obser- 
vations of galaxy clustering, have helped in turning cosmology into a precision science. The combined set 
of measurements point to a flat, low density (Cl m c± 0.3), accelerating Universe dominated by dark energy 
(Ode — 0.7). Modern redshift surveys aim to complement the above information about background cos- 
mological parameters by providing a detailed map of the distribution of baryonic matter on large scales. 
Redshift surveys have so far been primarily used for two purposes: (1) to get independent constraints on the 
matter density and several other cosmological parameters (see below) and (2) to measure the abundance and 
clustering properties of a variety of galaxies. The former exercise provides important consistency tests for 
cosmic parameters estimated using other methods, whereas the latter gives precious feedback to the physics 
of galaxy formation. In addition, it is also important to test whether the observed large-scale clustering 
properties of galaxies are consistent with the predictions of the preferred cosmological model(s). In order 
to carry out such a test, one needs to quantify the information content hidden in redshift space distribu- 
tion of galaxies. In this article we review the methods developed to study the geometry, connectedness 
and morphology of large scale structure (LSS) in redshift surveys, which can eventually help us meet such 
a goal. Such an approach gives valuable insights into the large-scale distribution of matter, and provides 
a framework within which to test the paradigms of gravitational instability and of the Gaussianity of the 
primordial density field. 

Impressive redshift surveys have been undertaken recently to map out the 3— dimensional distribution 
of galaxies. In the next section, we shall outline the history of redshift surveys, and focus on two major 
recent surveys, 2dFGRS and SDSS. LSS in redshift surveys exhibits phenomenally rich and complex texture, 
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wherein sheet-like and/or filamentary superclusters intersperse with large, empty regions called voids. Several 
salient aspects of the LSS of the Universe will be dwelt upon and we shall stress why it is important to study 
the morphology of LSS 2]. In subsequent sections we shall show that the complex morphological features 
of LSS are objectively quantifiable and that morphological information supplements traditional approach of 
understanding LSS using its n— point correlation functions and can be employed to confront theories with 
observations in an integrated manner. Statistical methods which probe the geometry and topology of LSS 
include percolation analysis 0, Q , genus curve 0, Q , minimal spanning trees 7] , void probability function 
8] , Voronoi tesselations 9] , as well as more recently introduced Minkowski Functionals (MFs) 0, ^| . For 
excellent reviews of these different approaches, the reader may refer to [2, EJ 0, • In this article, our 
focus will be on MFs. 



II. REDSHIFT SURVEYS 

Ever since sky maps of galaxy counts, e.g., Lick Catalogue revealed a rich pattern in LSS of galaxies 
in projection, it was realised that complementing this information with the redshifts of individual galaxies 
was essential to get a fully 3— dimensional view of the Universe. In such an exercise, one turns the redshift 
of a galaxy into its distance with the help of Hubble's law and is thus able to study the redshift space 
distribution of galaxies. 

Various small surveys were carried out in the late 1980's and mid 1990's in a controlled fashion; their 
positive scientific outcome providing impetus for the initiation of bigger ones. Of these perhaps the most 
influential was the Center for Astrophysics Survey pj|. This survey had an angular coverage of 6° x 120° 
and a depth of ~ 150/i _1 Mpc. The survey measured redshifts of about 2400 galaxies. The famous Great 
Wall centered around the Coma cluster was the most striking feature of this survey. The slice showed very 
clearly the "bubbly" nature of LSS with voids having sharp boundaries; the largest void having a diameter 
of ~ 50/i _1 Mpc. Subsequent surveys including follow-up CfA slices and the ESO Southern Sky SurvevfT^ 
amply confirmed the impression given by the first CfA slice. Indeed, the Southern Sky Redshift Survey 
(SSRS) discovered a second Great Wall in the southern sky. The emerging picture of filamentary and sheet- 
like superclusters |2(j encircling frothy voids prompted cosmologists to seek an explanation from theory, and 
this interplay between theory and observations is only now maturing into a truly quantifiable discipline. 
Another important survey, the Optical Redshift Survey (ORS)[2l| had a depth of 80/i _1 Mpc, but attempted 
a complete sky coverage (except for the zone of avoidance) . This survey measured 8500 redshifts in total, 
and was heavily used to estimate luminosity functions, galaxy correlations, velocity dispersions etc. In the 
mid 1990's, the Stromlo-APM redshift survey [2^ and the Durham/UKST redshift survey 2 3] were focussed 
towards the southern sky and led to fruitful results on correlation functions in real and redshift space, power 
spectra, redshift distortions, cosmological parameters, bias etc. However, the galaxy redshifts in all these 
surveys never exceeded ~ 10000. 

The Las Campanas Redshift Survey (LCRS) [2^ consisted of six 1.5° x 80° slices (3 each in the northern 
and southern galactic hemisphere), went to the depth of ^750 h~ 1 Mpc (z ~ 0.25) and recorded redshifts of 
about 25000 galaxies. LCRS presented an order-of-magnitude improvement over earlier surveys. This was 
the first deep survey of sufficient volume which could test ideas such as the statistical homogeneity of the 
Universe on large scales. Other important results from analysing LCRS included the luminosity function, 
second and third-order correlation functions, power spectra etc. LCRS is perhaps the only large survey 
wherein both the topology and morphology of LSS have been extensively investigated [25|, |2y, [23 and used 
to confront observations with theoretical models [2£| . 

A. The 2 degree Field Galaxy Redshift Survey (2dFGRS) 

The 2dF multi-fiber spectrograph on the 3.9m Anglo-Australian Telescope is a very impressive facility 
which can obtain spectra for up to 400 objects simultaneously over a field of view with 2° diameter. The 
spectrograph was put to use to sample two contiguous constant-declination strips of the sky, one in the 
northern galactic hemisphere [with angular coverage of 10° x 90° in (£, a)] and the other in the southern 
galactic hemisphere [with angular coverage ~ 15° x 90° in (£, a)]. Working in the Bj band, the survey reached 
a depth of z = 0.3 with median redshift of 0.11, and has compiled redshifts of about 2.5xl0 5 galaxies which 
are brighter than an extinction corrected magnitude of bj = 19.45 29]. The survey is already complete and the 
resulting galaxy correlation functions, redshift distortions and pairwise velocity dispersions [3(| demonstrate 
the superb quality of the data set. The science output of this survey and its overall cosmological relevance 
have been excellentlly reviewed in |3l| . We shall discuss here some of the results which are most relevant to 
our present discussion. 
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FIG. 1: Flux- limited 2dFGRS slices in northern and southern galactic hemisphere are shown. The supercluster- 
void network visually stands out, and can be properly understood by computing either a hierarchy of correlation 
functions or by alternate, geometric and topological quantifiers. In this article, we review methods to delineate 
individual superclusters and voids and describe their properties by using topological and geometrical diagnostics such 
as Minkowski Functionals and Shapefinders. 



One of the main science goals of 2dFGRS was to get improved estimates of the 3-dimensional power 
spectrum and the two-point galaxy-galaxy correlation function [32^. Since 2dFGRS sampled a cosmological 
volume as large as 0.1 Gpc] 3 , it became possible to measure these quantities over a large range of 
scales 0.01 < k < OAh Mpc -1 with an unprecedented accuracy [23, L13- Assuming a constant bias between 
dark matter and baryons on scales k < 0.15h Mpc -1 , Percival et al.(2001) compared the shape of the 
power spectrum with the primordial power spectrum predicted by linear perturbation theory of structure 
formation 36]. This allowed them to constrain the values of Q m , Qb after priors on the power spectrum index 
n and the Hubble constant h had been applied (Q m = 0.26 ± 0.03, Qb — 0.044 ± 0.016). Apart from this, 
2dFGRS was employed to obtain estimates for the luminosity function of early and late type galaxies [3^, 
and to study the luminosity dependent clustering of galaxies as measured using the correlation function 
statistic |22|. The latter study provided clues about the large scale bias between baryonic and dark matter 
and revealed that the former might follow the latter through a linear relationship (also see |4(l|). 

Cosmological parameter estimation is clearly one of the key results gleaned from large redshift surveys 
such as 2dFGRS. However, in our view several assumptions made in the process remain to be adequately 
justified, for instance, the assumption of linear, constant bias used in the exercise mentioned above is clearly 
inadequate. In fact, the study of biasing of galaxies with respect to the underlying mass distribution has 
its own history, and the issue has not yet settled completely. Below we shall briefly discuss the currently 
established ideas about bias. 

Figure ^ shows the distribution of galaxies in the flux-limited 2dFGRS slices. To the eye, galaxies 
are distributed anisotropically, with sharp qualitative similarities with the distribution of dark matter 
in N— body simulations (see below). The 3— dimensional power-spectrum, in a more appropriate form 
A 2 (k) = dP(k)/dln(k), measures the amount of power per unit interval in the logrithmic bin around 
wavenumber k. In simple terms, this is a statistical measure of the density fluctuations within spheres 
of radius R ~ 1/k. Clearly, in a measure such as this, the information on length-scales r < R is averaged out. 
Hence, one learns relatively little about the nature of the actual distribution of galaxies. Describing (and 
explaining) the visually rich filamentary pattern of superclusters of galaxies and quantifying the complex 
nature of the LSS is another, equally important application of redshift surveys. This problem has not been 
adequately addressed despite considerable recent progress 0, 0> El Ei 03 E3 • Clearly, the contribu- 
tion of higher order correlation functions is an important ingredient in such a program. Since higher order 
correlations are difficult to measure in practice, it would be useful if we could quantify the connectedness and 
non-Gaussianity of LSS as well as the shapes and sizes of individual superclusters and voids using geometrical 
diagnostics. The methods reviewed in the present article are intended to fill this gap in our knowledge and 
to complement more traditional, correlation function based approaches in developing a quantifiable picture 
of the LSS and its properties. 
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B. Sloan Digital Sky Survey (SDSS) 

The Sloan Digital Sky Survey (SDSS) is derived from a dedicated 2.5m telescope. The initial photometric 
program has been to measure the positions and brightness of about 10 8 objects in tt steradian (almost l/4 th ) 
of the sky. The survey is centered on the northern galactic pole and has an elliptical angular coverage of 
130° x 110°, where the semi- major axis runs along a line of constant right-ascension. 

Follow-up spectroscopy is planned to give redshifts of about 10 6 galaxies and 10 5 quasars. At 
the time of writing, the SDSS team has made its 2 nd data release, available for free download [48^ . 
It contains redshifts of ~ 3 x 10 5 galaxies, with the farthest galaxy at a redshift of ~0.3. (See 
http : //cas . sdss . org/astro/en/tools/search/SQS . asp.) 

The limiting magnitude of the SDSS survey is in the red band (ru m = 17.77), and its depth is similar 
to 2dFGRS. The survey will cover a cosmological volume as large as Gpc] 3 . Taking the length scale 

of 100 h~ 1 Mpc as a tentative scale of homogeneity, this is approximately thousand times a representative 
volume of the Universe. Hence when fully complete, SDSS can be fruitfully used to evaluate the cosmic 
variance for many statistics associated with LSS. The data are already large enough to permit accurate 
estimation of power spectrum 49], 3— point correlation function |Hoj . and of cosmological parameters [H"T| . 
Complementary studies of LSS in SDSS have been undertaken using MFs by |44| and using the minimal 
spanning tree- formalism bv 1491 . We shall mention later some results from an MF-based study of LSS in 
mock catalogues of SDSS |52 |. 

In order to make full use of the information encoded in the observed LSS, it is essential to understand the 
relationship between the number density of galaxies and the underlying mass density field. Since 1980s, it has 
been realized that the visible density field due to galaxies may be biased relative to the mass-density. In order 
to appreciate the implications that bias may have in the program of confronting theory with observations, 
below we briefly review the state of our knowledge in this regard. First we shall discuss the theoretical 
modelling of LSS using numerical simulations. 

III. N-BODY SIMULATIONS 

Strong evidence exists to suggest that dark matter is ~10 times more abundant than visible baryonic 
matter. Dark matter is usually assumed to be collisionless and to interact only through gravitation. Due 
to its relative abundance dark matter is clearly more important than baryons from a dynamical point of 
view. The best way to understand the evolution of a gravitating system consisting of N dark matter particles 
would be to actually solve for the force acting on a given particle due to the remaining particles. In a rather 
different form, by solving for the Poisson's equation at every stage of evolution of the matter-field, this is 
precisely the goal which modern versions of N— body simulation codes achieve. (For an excellent introduction 
to N— body simulations, see p53l IHEl 156^. 

Simulations allow us to probe the strictly nonlinear regime of the dark matter density field, and can provide 
us with detailed predictions of the dark side of the Universe at today's epoch. (The resulting dark matter 
distribution can be probed through observations of gravitational lensing - an area of great future promise.) 
The large scale structure distribution in dark matter could very well resemble that shown in Figure O We 
note from this figure, how the matter condenses along filamentary /sheet- like ridges between clusters and how 
voids get progressively more evacuated with the passage of time. Due to web-like and frothy appearance, 
LSS has been described at various times as being a Cosmic Web 57] or as Cosmic Froth [T^. 

During the past two decades the advent of parallel supercomputers have ensured a phenomenal improve- 
ment in the computing power available to tackle complex astrophysical processes. Hence, it has become 
possible to include complex gas-dynamical processes into simulations. Various recipes for cooling, star for- 
mation and feedback due to supernovae and galactic winds have been incorporated to mimic the formation 
and evolution of galaxies in a cosmological setup with simulated volumes ranging from 50 to 100 h~ 1 Mpc 
p58l [H^ . This has made it possible to develop mock catalogues of galaxies inspired by a given cosmological 
model. Relatively less time-consuming semianalytic methods have also been proposed which hope to achieve 
similar goals while bypassing the complex baryonic physics |l2l . 16(1 l6ll l62l |63| . Such mock catalogues can 
eventually be compared with clustering properties of galaxies derived from redshift surveys with the help of 
methods reviewed here. 



IV. THE NATURE OF BIAS 

Above we summarised the theoretical efforts being made to understand the clustering of dark matter. 
From the predictive point of view, given the fully evolved nonlinear mass-density field (both in the real 
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FIG. 2: The figure illustrates the development of nonlinear structures in a hydro dynamic N— body simulation from 
z — 6 (leftmost cube) to z — (rightmost cube) via z — 2 (middle cube). A near featureless density field evolves to 
produce filamentary and sheet-like superclusters which percolate through the box-volume. These are separated by 
large voids. Due to its web-like appearance, LSS is dubbed as the Cosmic Web. [Figure courtesy: Volker Springel; 
see http :/ / www. mpa-g arching, mpg. de / galform / data. vis / index, shtml] 



universe and in N-body simulations), we would like to know the exact locations where galaxies are formed. 
An important question related to this is whether the clustering properties of the resulting galaxy distribution 
(quantified in terms of n— point correlation functions or some other alternate statistics) are the same as the 
underlying matter-field, or whether the two are different. The answer to this, clearly depends upon the 
recipe which we utilise in grafting galaxies on the density field. 

During the mid-1990's it was shown that AGNs, IRAS-selected galaxies and galaxies selected in the optical 
bands cluster differently Earlier (during the 1980's) it was established that clusters of galaxies are more 
strongly clustered than galaxies themselves [g^. Recently it was proven that the more luminous galaxies 
are more strongly clustered [23,1^. Since these different visible tracers of LSS appear to cluster differently, 
their distributions are clearly biased relative to each other. There is hence, a reason to believe that the 
visible matter as a whole will be biased relative to the underlying dark matter. This phenomenon is termed 
as bias. Knowledge of the nature of bias will clearly enable us to connect the two distributions — visible and 
dark; observed and simulated. 

In cosmology, it is a conventional practice to quantify the clustering of LSS using the two-point correlation 
function or, equivalently, the power spectrum. Both these statistics are known to have a power-law form 
for galaxies |3E3- However, the correlation function of dark matter was shown to differ from ^ gg (r) in a 
complicated, scale-dependent fashion for various CDM-based cosmogonies 68, KM l70j . This suggests that the 
bias cannot be scale- independent, i.e., light may not follow the mass in a linear fashion. Recent investigation 
involving the three— point correlation evaluated for the SDSS dataset implies that galaxy biasing could be 
quite a complex and nonlinear process; see 50]. 

Kaiser (1984) explained the strong clustering of clusters using the high peak model [23. This model states 
that a rare, high-density fluctuation corresponding to a massive object of a given size, collapses sooner if 
it lies in a region of a larger-scale overdensity. This is because, the small-scale overdensity is aided by the 
surrounding large-scale overdensity so that the collapse-criteria is fulfilled relatively easily for such a halo. 
This "helping hand" from the long-wavelength mode means that overdense regions contain an enhanced 
abundance of massive objects with respect to the mean. Later this model was generalised to incorporate 
objects of any mass [7lll72( . 

The high-peak model is based on kinematic premises. It traces back the observed clustering of various 
classes of objects to the in-built structure in the primordial de nsity field. It could fruitfully be applied to 
study the clustering of galaxies only at high redshift (z ~ 2.5 — 3) |63,lzillz3, for the galaxy-scale fluctuations 
were prone to collapse at such early times. At such high redshifts, galaxies were shown to be highly biased 
relative to background 73]. Starting from z> 3, until the present, galaxies have very likely collapsed in all 
environments. The u nderlying ma ss-fluctuations have also grown. Thus, one predicts that the large bias 
(defined as h(r) = \/ igir) / ^dmir)) which existed earlier, should asymptotically approach unity on large 
scales. Lahav et al.(2002) showed for 2dFGRS data that the bias averaged over all scales is statistically 
consistent with a scale-independent, large-scale bias of order unity. However, we should see whether this 
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knowledge helps us achieve our original goal of connecting the simulations with the observed Universe, i.e., 
in populating the mass-density field with galaxies in the correct manner. 

A recent, Halo Model based approach is a fruitful step in this direction. Here all the complications of 
galaxy formation are encoded via the halo occupation number: the number of galaxies found above some 
luminosity threshold in a virialised halo of a given mass . In a nut-shell, the halo model describes nonlinear 
structures as virialized halos of different mass, placing them in space according to the linear, large-scale 
density field [t3 | . This model is analytically tractable, and correctly incorporates the dependence of bias on 
the 2 nd order correlations in the density field. Consequently, it correctly reproduces the observed 2— point 
correlation function and the power spectrum. Figure |3| compares the exact non-linear mass-density field 




FIG. 3: Left panel: Simulated dark matter distribution which shows complex morphology. Right panel: Halo Model 
corresponding to density field in the left panel, which successfully reproduces galaxy-galaxy correlation function and 
the power spectrum. The two are manifestly different, which points out the role which higher order correlation 
functions play in defining the cosmic structure. Figure courtesy 74]. 



with the halo model representation. The halo model is remarkably successful in obtaining correct correlation 
function, power spectrum etc. However, as shown in Figure 01 it is manifestly different from the fully evolved 
density field. This brings out an important aspect relevant to further study of LSS: the low-order correlations 
cannot be relied upon to provide us with a complete picture of clustering. In addition to the large halos of 
clusters and groups, today the rich filamentary /sheet like patterns in matter density field are also populated 
by galaxies [73, and one must theoretically understand the emergence of halos along filaments/sheets in order 
to complete the picture. In fact, the following experiment can be carried out: the galaxies can be selected in 
z = mass-density fields from two rival cosmogonies, so that the observed galaxy-galaxy correlation function 
£>gg( r ) an d the power spectrum are reproduced. In comparing our predictions with the LSS-data, if we 
confine ourselves merely to two-point correlation function, such an algorithm would validate both models to 
be good representatives of our Universe. Evidently this is not true. Clustering properties of galaxies in the 
two models would of course be different from each other (and perhaps even from LSS in the Universe). This 
serves to establish that the large-scale distribution of galaxies cannot solely be described by the two-point 
correlation function alone 76] . The resulting galaxy distributions can be shown to differ from each other by 
employing statistics which are sensitive to higher order correlation functions. Following this, we arrive at 
an important conclusion: the two distributions of galaxies must be compared with the help of higher order 
correlation functions or statistics derived therefrom. 

To summarise, although the high peak model ^ gives us promising results for galaxies assembling at high 
redshift and for Abell clusters at present, the predictive power of halo model in populating mass density field 
with galaxies is clearly inadequate (in fact, due to a variety of reasons, the bias could even be stochastic; see 
[73 )• In particular, for a combination of a cosmological model and a biasing scheme to be proven correct, 
it is not sufficient that the two explain the two-point correlation function for galaxies alone. In order to 
effectively compare the theoretical predictions for galaxy-distribution with the data from redshift surveys, 
one should include information from higher order correlation functions, and/or alternate statistics which can 
effectively quantify the complex morphology of the galaxy distribution. One such class of diagnostics the 
Minkowski Functionals - are introduced in the next section. 
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V. 



MINKOWSKI FUNCTIONALS 



A. Minkowski Functionals in Cosmology 



It is well known that, in contrast to a Gaussian random field (hereafter, GRF), a fully evolved nonlinear 
density field cannot be fully quantified in terms of its two-point correlation function; the latter being sim- 
ply the lowest and first of an infinite hierarchy of correlation functions describing the galaxy distribution. 
Furthermore, the bias in galaxies viz a viz dark matter can be nontrivial in nature, and can give rise to 
different clustering properties of galaxies and the underlying mass density field. Consequently the correla- 
tion function £ >gg (r) determined either for galaxies or for dark matter does not on its own validate a given 
cosmological model unless support is provided by other, statistically indepdend measures of clustering. It 
therefore becomes a challenge to compare the theoretical predictions for LSS with real observational data. 

If we knew all n— point correlation functions, we would have a complete description of the galaxy clustering 
process. However, estimating £ gg (r) for a sample of N galaxies requires knowing all pairs of galaxies in this 
sample, whereas calculating the 3-point function implies taking all triplets. The amount of computation 
escalates rapidly with n - the order in the n— point correlation function - especially for a large value of 
N. Besides, it is difficult to extract intuitively useful information from these statistics. The present article 
focuses on another class of statistics which complement the correlation function approach and which have 
the advantage of providing a physically appealing interpretation for the evolving density field. These are the 
Minkowski Functionals. 

For an excursion set involving particles embedded in n— dimensions, the Minkowski Functionals (hereafter, 
MFs) are defined on an (n — l)-dimensional hypersurface. There are n + 1 MFs in n dimensions. In this 
article, we shall be concerned with the 3-dimensional distribution of dark matter and/or galaxies. Hence, 
MFs will be defined on a 2-dimensional surface and they will reflect the physical properties of this surface 
(in the given instance , an isodensity contour referring to a supercluster of galaxies or a void). For a given 
surface the four MFs [HI EE Ell are, respectively 

1. Volume V, 

2. Surface area S, 

3. Integrated Mean Curvature C, 



Ri and R2 are the two principal radii of curvature of the surface in a given local neighbourhood and 
the integral is taken over the entire closed surface. 

4. Integrated Gaussian Curvature (or Euler characteristic) x> 



A related quantity which is more popular in Cosmology is the genus G. The genus is related to x by 



The 3-dimensional genus of an object is a topological invariant. It can be interpreted in terms of the 
connectivity of the surface. In simple terms it can be viewed as the number of independent cuts which 
one can make to the surface without breaking it into two separate pieces. A torus can be cut once and 
yet remain in one piece. Hence, its genus is 1. However, a sphere would separate into two pieces if a cut 
were made. Hence, its genus is 0. Two objects with the same value of genus are topologically similar: 
one can be obtained by continuously deforming the other. Thus, a sphere and a cube are topologically 
equivalent. More concretely, genus of an object is the number of handles that the object has, in excess 
of the number of holes which it encloses (e.g., see |8Cj). Thus, 




(i) 




(2) 



= l-±. 
2 



(3) 



G = [# of handles to the surface] — of holes enclosed by the surface]. 



According to this definition, a sphere has no handle, a torus has a single handle (equivalent to a sphere 
with one handle), and a pretzel has two handles. Introducing a hole or a bubble inside the surface 
reduces its genus by 1, whereas adding a handle to the surface increases its genus by 1. 
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An attractive feature of MFs is that they depend upon the entire hierarchy of correlation functions • 
Thus, MFs can indirectly help us capture the effect of n— point functions. Higher order n— point functions 
gradually become important as a primordial, featureless GRF evolves to develop nonlinear structures. This 
reflects in the behaviour of MFs for the corresponding density field. MFs are additive in nature, i.e., they 
can be studied for individual objects (say, clusters or superclusters defined using some prescription) as well 
as for the entire ensemble of such objects. 

Of the 4 MFs listed above, the genus G was already known to the cosmology community 0, E1E2- Tools 
like percolation analysis which can also be related to MFs, had earlier been introduced in cosmology by [8^ . 
There were strong reasons why this was the case. 

It is well known that a system evolving under gravitational instability becomes progressively more non- 
Gaussian. (The hypothesis of a primordial spectrum of density perturbations which is distributed in the 
manner of a Gaussian random field is supported by observations of the cosmic microwave background carried 
out by WMAP and other experiments.) For CDM-like cosmological models this results in more matter being 
concentrated in filamentary and pancake-like distributions whose coherence scale evolves with time. As a 
result the filling fraction at percolation progressively decreases, from 16% for the initial Gaussian random 
field to - 5% for LCDM |H El El H3 This implies that more matter is being transferred to regions which 
occupy a progressively smaller amount of space and are also spatially anisotropic and therefore percolate 
easily. Figuratively this corresponds to the prominance acquired by superclusters of galaxies over time - 
see Figure This theoretical picture which described the emergence of structure in an initially featureless 
medium agreed well with what observations were telling us about our universe. The emerging picture of LSS 
taken from redshift surveys shows us that galaxies are distributed preferentially along filaments and sheets 
which seem to encircle vast voids. The distribution of galaxies does not occupy much volume, and because 
of the presence of superclusters, percolates easily. Percolation theory can therefore be used to describe an 
important quality of the observed galaxy distribution - its connectivity. At the percolation threshold one 
can visit the far corners of the cosmic web by following filamentary /sheet-like over densities. The cosmic web 
is therefore a connected structure and this is clearly revealed in surveys such as Figure O Connectivity is of 
course a mathematical notion, and it is possible to precisely quantify connectivity of an object or a system by 
evaluating its percolation properties and its genus. Doroshkevich(1970) gave analytic formula for the genus 
of a GRF which showed that at the median density threshold, GRF would exhibit a sponge-like topology 
[82^ . Simulations of LSS within CDM-cosmogonies provided similar visual impression when smoothed on 
sufficiently large scales. Like percolation, the genus too has come to be considered as a useful probe of 
the initial Gaussianity (as well as the final non-Gaussianity) of the cosmic density field; see for instance 

An important indicator of non-Gaussianity in a distribution is the percolation curve. The percolation curve 
describes the volume fraction (or filling factor - FF) in the largest structure (cluster /void) as a function of 
the density contrast threshold 5 [&|. Salient features of the percolation curve are illustrated in Fig. 0]in 
which density fields (with an initial power spectrum P(k) ~ k~ 2 ) evolve from an epoch when the scale of 
nonlinearity is /cnl = 64fc/, 16fc/, 4fc/, (kf is the fundamental mode corresponding to the box-size of our 
N-body simulation). In Fig. 0] percolation curves for clusters (thick solid lines) and voids (thick dashed 
lines) are shown as functions of the density contrast 5. Starting from a high density threshold (small FF) we 
find several isolated clusters (corresponding to peaks of the density field). Lowering the density threshold 
further results in the merger of clusters leading to a rapid growth in the percolation curve and to the onset 
of percolation. A further lowering of the threshold to very small values results in the merger of almost all 
clusters so that FF — > 1. An identical procedure followed for underdense regions by gradually increasing 
the density contrast threshold (increasing FF) results in a similar functional form for the volume fraction 
in the largest void. In our samples the largest cluster percolates between opposite faces of the cube when its 
filling factor is about half the total FF. In most cases percolation also coincides with the highest jump in 
the volume of the largest cluster which was used by 89] as a working definition of the percolation threshold. 

The solid/dashed vertical line in Fig. ^represents the density contrast threshold 5c below/above which 
clusters/ voids percolate. We note that both clusters and voids percolate over a range of overlapping density 
contrasts — a feature that is only possible in three or more dimensions — and corresponds to what is 
commonly called a sponge topology for the density distribution. 

As the simulation evolves 5 c increases monotonically for n — —2, as structures form and align on in- 
creasingly larger scales. From Fig. 0]we see that voids find it easier to percolate as the simulation evolves, 
as a result the range in densities when both phases percolate initially increases, enhancing the extent of 
sponge-like topology in the distribution. 

The right panel in Fig. 21 shows the evolution of the genus curve for the same simulation. The top panel 
corresponding to an early epoch shows the genus curve largely retaining its bell-shaped form which it has 
for the primordial Gaussian random field G{y) = Ail — v 2 ) exp(— ^ 2 /2) (IHESl). During later epochs the 
density distribution grows progressively more non- Gaussian and this is reflected in the change in both shape 
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FIG. 4: Percolation (left panels) and genus (right panels) curves are plotted as functions of the density contrast S for 
scale free models of gravitational clustering with the perturbation spectral index n — — 2. Solid and dashed curves in 
the left panels correspond to the percolation curve for the largest cluster and void respectively. Vertical solid/dashed 
lines mark the threshold describing percolation between opposite faces of the cube for clusters/voids respectively. 
Figure courtesy of |85|. 



as well as amplitude of the genus curve in Fig. 21 

Although a study of the genus gives us a useful handle on the connectivity of a density field, we still lack 
information about the morphology of LSS. For example, if the domain of information is solely restricted to 
knowing genus, a filament with one handle would be considered to be identical to a sphere with one handle 
or a pancake with one handle. Therefore, to gain objective insights into the nature of the supercluster-void 
network, we must complement genus of an object with quantities which in some way, characterise geometry of 
that object. This is precisely the role played by the first three MFs, Volume V, Surface Area S and Integrated 
Mean Curvature C. These MFs change when local deformations are applied to the surface, and hence, these 
may be useful to glean information about the typical size of an object. The final goal here would be to use 
this geometric information and give an objective meaning to generic sizes and shapes of super clusters and 
voids belonging to LSS. In this article we will show that MFs can be utilised to characterise the geometry and 
topology of the cosmic density field, and eventually can be employed to get information about the morphology 
of the supercluster-void network in LSS. Because MFs depend on the full hierarchy of correlation functions, 
we can use them to compare and distinguish two rival cosmological models. Recently, Matsubara (2003) 
derived semianalytical expressions of MFs for weakly nonlinear cosmic density field by using a second-order 
perturbative formalism |80(. The results obtained by Matsubara can be used in conjunction with the methods 
reviewed here to obtain insight into the dynamics of the clustering process of a density field smoothed on 
sufficiently large scales. The role of bias on large scales can also be directly probed using such methods. 



B. Morphology of LSS with Shapefinders 

Sahni, Sathyaprakash & Shandarin (1998) showed that the MFs could be used to evaluate both the size 
as well as the shape of a three-dimensional object such as a supercluster or a void l?)^. The size is given in 
terms of three Shapefinders which have dimensions of length. These are defined as ratios of MFs, and are 
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conviniently termed as Length (£), Breadth (B) and Thickness (T). 

£ = < 6 » 

The three Shapefinders are defined using spherical normalization, so that for a sphere of radius R, T = 
B = C = R. 

The above measures quantify the size of the object in question. To further quantify the shape of these 
objects, two dimensionless Shapefinders have been defined as follows. These are Planarity (V) and Filamen- 
tarity (J 7 ). 

B — T 

V = - — - (7) 

s + r 1 ; 

T = cTb- (8) 

Sahni et al.(1998) showed that (V^J 7 ) quantify the shape of both simple as well as complicated objects. 
Thus, for an oblate ellipsoid its planarity is relatively large and its filamentarity is small. For a prolate 
spheroid, the reverse is true. The statistics respond monotonically to deformations of these surfaces. 

If we include the genus then, the triplet of numbers (V^F, (?), can define a three-dimensional Shape- 
Space which may be used to represent the distribution of shapes of superclusters and their topologies in a 
given cosmological density field. Note that the inclusion of genus is significant, since most superclusters at 
moderate density thresholds can have a spongy texture with several loops and/or branches emanating from 
a central body (see Figure fTTI for an illustration). 



C. Methods of Estimation 



In cosmology the distribution of galaxies in space can frequently be regarded as being an example of a point 
process. There have been several attempts of evaluating MFs for a given point process. For example, the 
Boolean grain model [lQ, ti£j| would decorate the input set of points with spheres of varying radii, and would 
study the morphology of the structures which contain overlapping spheres. According to another approach 
[Tl| one smooths the distribution of points using a suitable window function, and defines the density field 
on a grid. MFs can be evaluated for connected structures defined on a grid by using Crofton's formulae or 
Koenderink Invariants [Tlj . 

A new approach for determining MFs has been reported in 93]. Rather than working with a grid- 
approximation for connected structures, these workers go a step forward, and model surfaces for these 
objects. Minkowski Functionals are evaluated for the isodensity surfaces online, i.e., while these are being 
modelled using an elaborate surface triangulation technique. Sheth et al.(2003) describe the algorithms for 
surface construction and for the evaluation of MFs on these surfaces. Sheth (2004) discusses the method of 
implementation of these algorithms which has culminated into a robust and accurate software SURFGEN 
(short for "SURFace GENerator" ) 94] . SURFGEN has been tested for accuracy against Gaussian random 
fields and against a variety of simply and multiply connected eikonal surfaces. SURFGEN evaluates the 
geometric MFs to better than 1% accuracy in most cases, while the genus is evaluated exactly [93|. SURF- 
GEN has since been employed to study geometry, topology and morphology of LSS both in dark matter and 
galaxies within a variety of cosmogonies. (SURFGEN could also be relevant in areas of science where the 
properties of surfaces becomes important such as in medical imaging (tomography) and in condensed matter 
physics.) 

In the next section we shall discuss a set of new results obtained by applying MF-based techniques to 
LSS-data. 



VI. RESULTS 



A. Morphology of Large Scale Structure using Shapefinders 



Given a cosmic density field, Minkowski Functionals (MFs) can be evaluated for a supercluster (void) 
defined as a connected overdense (underdense) region above (below) certain physically motivated threshold 
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FIG. 5: Three global MFs for ACDM (solid lines), rCDM (dot-dashed lines) and SCDM (dashed lines) are shown as 
functions of the overdense volume fraction FFv- The three models have appreciably different morphology and hence 
can be distinguished from one another on the basis of MFs. For futher discussion, please refer to the text. Figure 
courtesy of |93j . 



of density. SURFGEN-like technique would isolate such a structure, and would model its isodensity contour 
in order to precisely estimate its MFs. MFs referring to a given structure are called partial MFs. By knowing 
partial MFs, one can further measure shape and size of the concerned structure using Shapefinders. Further, 
MFs are additive in nature. Hence, at any given threshold of density, global MFs pertaining to the entire 
density field can be obtained by summing over partial MFs of the constituent overdense/underdense objects. 
While partial MFs are useful to measure shapes and sizes of individual superclusters and voids, global 
MFs can be employed (1) to discriminate between models, (2) to probe the effect of nonlinear gravitational 
evolution on the density field and (3) to capture the nontrivial effects of bias on the structure of the cosmic 
web. We shall review here some of the promising results obtained in these directions. 



B. Discriminating between rival models of LSS 

The growth-rate of density perturbations differs between various cosmogonies. As a result, the clustering 
properties of the matter-distribution also turn out to be different. A serius candidate statistic of LSS should 
be able to capture such differences and discriminate between rival cosmological models. 

The three upper panels of Figure show clustered distribution of dark matter at z = in three rival 
cosmological models — SCDM, rCDM and ACDM respectively, simulated by Virgo group |z3|- The three 
models are simulated such that the observed abundance of rich clusters of galaxies can be reproduced in all of 
them. SURFGEN was employed to quantify the visually apparent differences between these models. Global 
MFs were evaluated at a set of density levels, and were studied with respect to FFv, the fractional overdense 
volume The lower 3 panels of Figure show variations in 3 global MFs — area, curvature and genus — 
as density level is progressively brought down. Notice that MFs for SCDM (dashed lines) show the largest 
amplitudes, those for ACDM (solid lines) show the smallest amplitudes with rCDM curves (dot-dashed lines) 
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falling between the two. Since the number of rich clusters and their relative spatial distribution in all the 
models is the same, this difference can be attributed to difference in matter-distribution along sheets and/or 
filaments connecting the clusters |93l |94[ . Global MFs are thus noted to be remarkable at discriminating 
between models (also see |95j). 

In Section |^ we noted that the distribution of galaxies cannot be fully quantified using its 2— point 
correlation function £ gg {r) alone. However, the galaxy distribution predicted by a given cosmological model 
must at least reproduce the observed € gg (r). Cole et al.(1998) generated mock catalogues of galaxies in a 
variety of cosmogonies subject to the above constraint. Thus, galaxies were selected from the simulated 
z = mass-distributions so that the observed Cgg( r ) 1S reproduced. The catalogues were further used to 
construct mock 2dFGRS and SDSS redshift surveys in such models [96|. Since MFs depend on the hierarchy 
of correlation functions, these can be relied upon in distinguishing the galaxy distributions resulting from 
different models. Sheth (2004) made such a comparison between mock SDSS catalogues of galaxies from 
ACDM and rCDM models by evaluating global MFs of density fields smoothed with a Gaussian window 6 
h~ x Mpc wide. Figure shows the global mean curvature and genus for the two models evaluated using 

I i i t | i i i | i I i | i i i | i i i | I I | i i i | i i i | i I i | i i i | i i i | 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

Volume filling fraction FF Volume filling fraction FF 



FIG. 6: Two global MFs — curvature and genus — are evaluated at 50 density levels at a common set of volume filling 
fractions and are studied here with respect to FFy. The values are normalised to the volume of [100 h~ 1 Mpc] 3 . The 
global MFs for rCDM model are averaged over 10 realizations (solid lines) and the errorbars represent la deviation. 
Assuming the same level of accuracy for ACDM, we may conclude that these MFs with volume parametrisation can 
indeed clearly distinguish rCDM from ACDM. Figure courtesy of |52|. 



SURFGEN. Note that the MF-amplitudes in ACDM are larger than those for rCDM model, and that the two 
models are successfully distinguished from each other using MFs 52], even in the case when the distributions 
share the same £ gg (r). (The rCDM curves are averaged over 10 realizations made available by |9q.) MFs 
can thus be reliably adopted to compare mock catalogues with the observed LSS in redshift surveys. 

C. MFs and the scale-dependent bias 

We reviewed the status of our knowledge about bias in Section HVl We raised an important question as 
to whether the clustering properties of dark matter and the biased galaxy distribution could be different or 
otherwise. This question was investigated by [52| by employing the mock SDSS catalogues of galaxies in 
ACDM and rCDM models. The clustering properties of galaxies were measured in terms of MFs. 

As noted by 0,113, the amplitude of the genus-curve drops as the N— body system develops phase corre- 
lations. Given two density fields, the system with larger genus-amplitude shows many more tunnels/ voids 
which are, therefore, smaller in size. As time progresses, the voids are expected to expand and merge, leading 
to a drop in the genus- amplitude, while the phase correlations continue to grow. With this simple model, one 
could correlate the amount of clustering with the relative smallness of the amplitude of genus, and therefore, 
of the area and the mean curvature. Going by this reasoning, we may infer from Figure that the ACDM 
galaxy-distribution is relatively more porous and would show less coherence on large scales compared to 
that due to rCDM. As illustrated in Figure [Sfsee |93|)- the dark matter distribution of rCDM model due to 
Virgo group shows considerably larger amplitudes for the MFs compared to the ACDM model, whereas, we 
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FIG. 7: Left panel: Global genus curves of ACDM and rCDM matter distributions are shown. The values refer to a 
volume of [100 /i _1 Mpc] 3 . The la error-bars are due to 5 realizations of both the models, each with 25% of the total 
number of particles. Calculations are carried out on a grid of resolution 3.5 h~ 1 Mpc after smoothing the density 
fields with 6 /i _1 Mpc Gaussian kernel. Right panel:The galaxy catalogues due to the two models are analysed in real 
space under identical conditions. This figure illustrates the phenomenon of phase-reversal. For more details, please 
refer to the text. Figure courtesy of |5^. 



find the reverse trend in the MFs of the galaxy distributions due to the same two models. Evidently biasing 
appears to be a source of this effect. 

To establish this effect more firmly, the rCDM and ACDM dark matter Virgo simulations were analysed 
by adopting the same resolution (£ g =3.5h~ 1 Mpc) and the smoothing scale (A s =6/i _1 Mpc) as utilised in 
the analysis of mock catalogues. So as not to introduce any bias due to redshift space distortions, global 
MFs of the galaxy catalogues were computed in real space. The effect of biasing is most dramatically seen 
in the global genus-curve of the dark matter and galaxies (see Figure |7J). We notice that the amplitudes 
Gdm < ^Lm> wnereas Qg > &g- Here subscripts DM and G stand for "dark matter" and "galaxies", 
respectively. We may conclude from here that the matter distribution in ACDM is less porous than that in 
rCDM, but this trend reverses as we investigate biased distributions of galaxies due to the same two models: 
the galaxy distribution due to ACDM appears to be more porous than that due to rCDM. This could perhaps 
mean that a simple scheme of scale-dependent biasing leads to different degree of phase-correlations in the 
galaxy distribution as compared to the underlying mass distribution. This can be dubbed as a phenomenon 
of phase- mismatch between the two distributions —visible and dark 52]. This is supported by the fact that 
|9q require large anti-bias in high density regions of the ACDM model. 

To conclude, the study of the global MFs reveals that, local, density-dependent bias could lead to an 
apparent phase-mismatch or different phase-correlations among the dark matter and the galaxies. It therefore 
becomes important to include realistic treatment of bias before comparing theoretical predictions about LSS 
with redshift surveys. We noted earlier that morphology of LSS appears web-like with frothy voids separating 
sheet-like and/or filamentary superclusters. The rich texture of LSS is a culmination of two interconnected 
factors: (1) strong non-Gaussianity induced in cosmic density field due to gravitational dynamics and (2) 
biased formation of galaxies relative to the underlying mass. Superclusters are thus strongly non-Gaussian, 
baryonic structures. Unlike galaxies and their clusters, superclusters are as yet unrelaxed and dynamically 
evolving; their morphological properties may be expected to be sensitive to the cosmological parameters 
of the Universe. It is hence important to quantify their shapes and sizes and relate them with physically 
relevant other quantities such as the mass enclosed or volume occupied by them. Similarly, voids too can be 
studied, and the combined properties of the supercluster-void network could provide us with yet another, 
independent check on the cosmological parameters. Since modern redshift surveys are deep and large enough, 
this may perhaps be the most effective method to test theoretical predictions. 

Recently, a comprehensive morphological study of the supercluster-void network in ACDM cosmogony was 
carried out by j87|. Dark matter distribution at z = was smoothed with 5 h~ x Mpc Gaussian window 
to produce a cosmic density field sampled over a box of size 239.5 IrT 1 Mpc. Superclusters and voids were 
defined as over dense and underdense connected regions, respectively and their MFs and Shapefinders were 
measured using SURFGEN. 
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TABLE I: The ten most voluminous superclusters (determined at the percolaton threshold) are listed with their 
associated Minkowski functionals (Volume, area, curvature and genus) and Shapefinders T, B, C in case of the ACDM 
model simulated by Virgo group. The first row describes the percolating supercluster and appears in boldface. It 
should be noted that the interpretation of C as the 'linear length' of a supercluster can be misleading for the case of 
superclusters having a large genus. In this case C X (G + 1) provides a more realistic estimate of supercluster length 
since it allows for its numerous twists and turns. 



Model Volume Area Curvature Genus Shapefinders (h x Mpc) 
(h^Mpc) 3 (h^Mpc) 2 h^Mpc T B £ 



ACDM 8.45 x 10 4 


4.5 x 10 4 


6.16 x 10 3 


6 


5.63 


7.30 


70.03 


5 = 2.31 3.21 x 10 4 


1.7 x 10 4 


2.33 x 10 3 


1 


5.63 


7.34 


92.74 


2.22 x 10 4 


1.22 x 10 4 


1.76 x 10 3 


1 


5.47 


6.93 


70.00 


2.13 x 10 4 


1.02 x 10 4 


1.29 x 10 3 





6.27 


7.90 


102.7 


1.21 x 10 4 


6.35 x 10 3 


8.78 x 10 2 





5.71 


7.22 


69.95 


1.19 x 10 4 


6.82 x 10 3 


1.02 x 10 3 


1 


5.25 


6.70 


40.53 


1.08 x 10 4 


6.12 x 10 3 


9.16 x 10 2 





5.3 


6.68 


72.90 


1.01 x 10 4 


5.46 x 10 3 


8.19 x 10 2 





5.57 


6.65 


65.24 


8.6 x 10 3 


4.74 x 10 3 


7.07 x 10 2 





5.44 


7.00 


56.33 


8.1 x 10 3 


4.25 x 10 3 


5.71 x 10 2 


1 


5.72 


7.44 


22.74 



Superclusters in ACDM model were found to percolate through the Universe while occupying fractional 
overdense volume as small as 7% 87]. This points to the connectedness of the overdense regions in the 
Universe. From their tendency to effectively percolate the volume, it is apparent that superclusters must 
show strong departure from sphericity. Voids, on the other hand, occupy as large a fractional volume as 
22% before the underdense regions percolate |87j . Evidently, the morphology of voids would be markedly 
different from superclusters. In shape, they are likely to be more isotropic. 

Shandarin, Sheth and Sahni (2004) further supported above inferences by studying the sizes of superclusters 
and voids with respect to the enclosed mass and occupied volume, respectively. Figure |HI shows the results. 
Notice that all three sizes of the structures show significant correlation with mass enclosed (for superclusters) 
and volume occupied (for voids): the larger the mass, the larger the size of SC and more the volume, the 
larger is the size of the void, as expected. The solid lines show the radius of a sphere having the same volume 
as a given object (R = (3F/47r) 1//3 ). The thickness and breadth approximately double their values and length 
grows by over an order of magnitude when the mass increases from about 10 14,5 M to 10 16,5 M@. Both 
the thickness and breadth are considerably smaller than the radius R of a sphere having similar volume for 
large superclusters (M> 10 15 Mq). On the other hand the length is considerably greater than R. This is a 
clear manifestation of the anisotropy of large-scale superclusters (see Figure [HJ. Voids demonstrate similar 
correlations with their volumes as shown in the right panels of Figure [SI Voids, especially large ones, are 
also anisotropic, e.g., see [87| for illustrations. It was shown for example, that both superclusters and voids 
show significantly larger filamentarity in proportion to the mass enclosed and volume occupied, respectively. 
Superclusters show negligible planarity whereas, voids could be relatively more planar |87j . 

In their visual appearance, filaments stand out remarkably well in our impression of LSS. However, the 
occasional occurrence of vast, sheet-like superclusters such as the SDSS Great Wall (which is almost twice 
as large as the CfA Great Wall ^JEll; see Figure EH) raises questions as to how significant are sheets with 
respect to filaments on large scales in the actual Universe. In future, if the sheet-like superclusters in redshift 
surveys are confirmed and their counterpart structures not found in N— body simulations, it may present a 
clear dichotomy between theory and observations. 

To investigate this question in the context of simulations, |93[ studied the correlation in the Shapefinders of 
ACDM superclusters identified at percolation. FigurelTUlis a scatter plot of Shapefinders T, S, C for clusters 
in ACDM defined at the percolation threshold. The strong correlation between T and B in the left panel 
indicates that two (of three) dimensions defining any given cluster assume similar values and are of the same 
order as the correlation length. (T~6~5/i _1 Mpc for the largest superclusters.) The clustering of objects 
near T~B(V — 0) in this panel suggests that the superclusters are either quasi-spherical or filamentary. 
The superclusters in simulations thus do not appear to be planar. The scatter plot between B & C in 
the right panel of Figure [ITU breaks the degeneracy between spheres and filaments. The mass-dependence 
of morphology is highlighted in this panel in which larger dots denote more massive objects. This figure 
clearly reveals that more massive clusters/superclusters are, as a rule, also more filamentary, while smaller, 
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FIG. 8: The length, breadth, and thickness versus mass for superclusters and versus volume for voids at percolation. 
Solid circles show the relation at percolation thresholds: FFc = 0.07 for superclusters and FFv = 0.22 for voids. 
Crosses show the parameters before percolation (FFc =0.06 for superclusters and FFv = 0.21 for voids) and empty 
triangles after percolation (FFc — 0.08 for superclusters and FFv = 0.23 for voids). Solid lines show the radius 
of the sphere having the same volume as the corresponding object. Note the logarithmic scale used for the length. 
Three lines correspond to three different thresholds. Figure courtesy of |87). 



less massive objects, are more nearly spherical 93]. 

For Shapefinder-like statistics, it is important that the measured properties of an object also conform 
with its visual appearance. Figure ITT1 shows the largest percolating supercluster in the ACDM cosmogony 
for which T = 0.81 |93|. Note that the supercluster is indeed highly filamentary in agreement with the 
measured, high value of T . 

Earlier we noted results from MF-based study of mock SDSS catalogues due to ACDM and rCDM mod- 
els. Global MFs were shown to successfully discriminate between the two models. Further, we noted the 
phenomenon of phase-mismatch between dark matter and galaxies due to biasing. In this subsection, we 
noted that superclusters in a typical CDM-cosmogony are generically filamentary. Could we use the length 
of the filamentary superclusters as discriminatory statistic between models? In other words, are the mor- 
phological properties of superclusters indeed sensitive to the underlying cosmic parameters? The 10 most 
voluminous superclusters in dark matter Virgo simulations of ACDM, SCDM and rCDM models were shown 
to have considerably different morphological properties. In fact, the percolating ACDM supercluster (shown 
in Figure ITTJl was found to be more filamentary and topologically more simpler (G = 6) than its counterpart 
structures in SCDM and rCDM (G ~ 20). 

To answer the same question in case of galaxy distributions, studied the abundance of superclusters 
larger than certain length C from within a given large ensemble of superclusters selected from mock SDSS 
catalogues of ACDM and rCDM models. This quantity is referred to as Cummulative Probability Function 
(CPF). Figure El shows the results, where the CPF(£) have been evaluated at FFy=7.8 per cent (at 
the onset of percolation) and at FFv =11. 7 per cent (after the percolation). Note that at the onset of 
percolation, the length of the longest rCDM superclusters could be as large as 90 /i _1 Mpc, whereas their 
ACDM counterpart structures, which exhibit same degree of statistical significance, are relatively shorter 
with £ max =55 h~ 1 Mpc. We conclude that the rCDM superclusters tend to be statistically longer than their 
ACDM counterpart structures. The large-scale coherence in the superclusters is attributed to the phase- 
correlations in the density field; the higher the degree of phase-correlations, the larger the length-scale of 
coherence. Based on this, and the results reported earlier, we would anticipate the rCDM superclusters to 
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FIG. 9: The Sloan Great Wall situated at z c± 0.08 is shown in the upper panel of the figure. For comparison, a 
relatively nearer structure of the CfA slice is also shown. The linear extent of the Sloan Great Wall is ~ 420/i _1 Mpc, 
which is about twice larger than the CfA great wall. Figure courtesy of 



be longer than those in ACDM. As we can see, the result reported here agrees well with this anticipation. 
This is a considerable success for the ansatz of the Shapefinder quantifying length, for as we noted above, it 
helps us capture the relative effect of phase-correlations among rival models of structure formation. 

D. Applications to Redshift Surveys 

We end this section by providing a brief summary of the results which have emerged out of the analysis 
of the observed LSS. We further outline the scope of these methods in light of the newly available large 
datasets. 

The methods of topological analysis of LSS have been well developed since late 1980s 6j. Hence, most 
of the investigations complementing the standard correlation function based approach are in the wake of 
studying the connectivity and topology of LSS. Almost all the topological analyses are carried out to test 
whether the galaxy-distribution is consistent with the underlying Gaussian, random-phase primordial density 
field. Topology of LSS revealed both in LCRS data 25] and SDSS data has been found to be consistent 
with Gaussian, random-phase hypothesis for the primordial density fluctuations. SDSS Early Data Release 
has been further found to agree well with the predictions of ACDM concordance cosmology model |4j1. 
Springel et al.(1998) reached similar conclusions based on their analysis of the 1.2 Jy redshift survey |97j . 
Connectivity of LSS of LCRS has been studied by 86]. These authors detected a network of filamentary 
superclusters in LCRS slices, which they found to be consistent with the gravitational instability paradigm 
of structure formation. 

Recently there has been an upthrust of interest in complementing earlier topological studies with geometry 
of LSS. Probably the first morphological analysis of LSS was done using 1.2 Jy redshift survey data by |9£j. 
MFs for SDSS early data release (EDR) have been evaluated by [44| who find the data to be consistent 
with the clustering predictions of the ACDM cosmogony |l00|. Basilakos (2003) has measured shapes of 
superclusters in SDSS (EDR) and finds them to be predominantly filamentary. Due to their large volume, 
SDSS and 2dFGRS can be fruitfully employed to test predictions for MFs in a weakly nonlinear regime. An 
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FIG. 10: Scatter plot for the pair of Shapefinders T , B (left panel) and B , C (right panel) defining the morphology 
of clusters/superclusters in the ACDM model. The strong correlation between T and B in the left panel near the 
line V — indicates that two of the three dimensions defining a cluster are equal and of the same order as the 
correlation length (p± few Mpc). Judging from the left panel we find that clusters/superclusters in ACDM are either 
quasi-spherical or filamentary (since both satisfy T~23=> V c± 0). The degeneracy between spheres and filaments is 
lifted by the right panel which is a mass-weighted scatter plot for the Shapefinders 23, C. Each dot in this panel refers 
to a cluster and its area is proportional to the fraction of mass in that cluster. The concentration of points near the 
line T — (B— C) reflects the fact that a large number of smaller clusters are quasi-spherical. The more massive 
structures, on the other hand, tend to be filamentary and the largest and most massive supercluster has T — 0.81. 
All objects are determined at the percolation threshold. Figure courtesy of |93|. 



interesting application has been illustrated by 1 1 ll | . who use mock SDSS catalogues (inspired by ACDM 
cosmology) for the purpose. 

Morphology of LSS in LCRS has been studied by [2(j who detect significant filamentarity in LCRS slices 
compared to an equivalent random distribution. The LSS-data from recent redshift surveys such as SDSS 
and 2dFGRS yet await an integrated morphological analysis and a detailed comparison with theoretical 
predictions. Meanwhile, encouraging efforts have been recently made to constrain the scale of homogeneity 
of the LSS by investigating the largest length-scale of filaments in LCRS slices |27j . Robust 2-dimensional 
morphological statistic such as the Shapefinder T have been used for the purpose. The scale of homogeneity 
is found to be c^70— 80 h~ x Mpc [2?| . ACDM mock LCRS catalogues have been further confronted with the 
LCRS slices using the same techniques. A uniform, moderate large scale bias of 1.15 seems to give good 
agreement between the model and the LCRS slices. 
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FIG. 11: The largest (percolating) supercluster in ACDM. This cluster is selected at the density threshold which 
marks the onset of percolation (S pe rc = 2.3). As demonstrated in the figure, the cluster at this threshold percolates 
through the entire length of the simulation box. It is important to note that the percolating supercluster occupies 
only a small fraction of the total volume and its volume fraction (filling factor) is only 0.6%. Our percolating 
supercluster is a multiply connected and highly filamentary object. Its visual appearance is accurately reflected 
in the value of the Shapefinder diagnostic assi| med to this supercluster: (T,#,£) = (5.63, 7.30, 70.03) /i _1 Mpc and 
(P,F,G) = (0.13,0.81,6). Figure courtesy of 
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FIG. 12: Shown here are the cumulative probability functions of length C for superclusters of galaxies in rCDM and 
ACDM at two thresholds of density corresponding to FFy=7.8 per cent (just before the onset of percolation) and 
FFv=11.4 per cent (after the onset of percolation). We find the CPFs of the two models to be distinctly different at 
longer length-scales. The dashed line refers to the CPFs due to the first realization of rCDM which shares the initial 
set of random numbers with ACDM. The largest superclusters of rCDM are systematically larger than those due to 
ACDM. Figure courtesy of [5^. 
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